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Abstract. Salmon's nearly geostrophic model for rotating shallow- water flow is derived in 
full spherical geometry. The model, which results upon constraining the velocity field to 
the height field in Hamilton's principle for rotating shallow-water dynamics, constitutes 
an important prototype of Hamiltonian balanced models. Instead of Salmon's original 
approach, which consists in taking variations of particle paths at fixed Lagrangian labels 
and time, Holm's approach is considered here, namely variations are taken on Lagrangian 
particle labels at fixed Eulerian positions and time. Unlike the classical quasigeostrophic 
model, Salmon's is found to be sensitive to the differences between geographic and 
geodesic coordinates. One consequence of this result is that the /3 plane approximation, 
which is included in Salmon's original derivation, is not consistent for this class of model. 
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1. Introduction 

The rotating shallow-water (SW) equations constitute a paradigm for geo- 
physical fluid motions ranging from fast timescale dynamics, associated 
with inertia-gravity waves, to slow advective-timescale dynamics, associated 
with nonlinear vortical motions and Rossby waves (cf. IGilll 119821 |Pedloskyl 
1987 ). This set of equations constitute the "primitive" equations on which 
different approximations are usually performed. In this paper I deal with 
those approximations which involve the introduction of balance relations or 
constraints that lead to filtering out the fast degrees of freedom. Terms com- 
monly used to denote the resulting models are "balanced," "constrained," 
or "intermediate;" the latter, in particular, reflects the fact of being at 
a level which is in between the primitive equations and the equations for 
geostrophic motion. For an extensive review on the wide variety of balanced 
models that exists in the literature the reader is referred to I Allen et al.l 
(199251. 



*To appear in O. U. Velasco-Fuentes et al. (eds.), Nonlinear Processes in Geophysical 
Fluid Dynamics, Kluwer Academic. 
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Of particular interest are those balanced models derived by performing 
approximations directly in Hamilton's principle (HP) for SW dynamics as 
proposed by Salmon(1983, hereafter referred to as S83). This procedure al- 
lows the fundamental symmetry-based conservation laws of the underlying 
primitive system to be preserved. The approach consists in substituting 
leading order balance relations and asymptotic expansions into HP before 
taking variations. In particular, S83's model is derived by constraining the 
velocity field to the height field in the form of a geostrophic balance relation, 
i.e. between the pressure gradient and the Coriolis force. This so-called 
LI model, however, was shown to produce less accurate solutions to the 
SW equations than those produced by other non-Hamiltonian intermediate 
models ( [Allen et al., 1990a| |Allen et al., 1990b| |Barth et al., 1990| >. This 
is indicative of the known fact that possession of Hamiltonian structure 
is no guarantee of model's accuracy. Nevertheless, other balance relation 
choices — potentially more accurate than that considered by S83 — are pos- 
sible JAlleri and Holm, 1996| |Allen et al., 2002| ). This fact makes the LI 
model an important prototype of constrained Hamiltonian models, and 
thus motivates the present study. 

The LI nearly geostrophic model, as well as its relatives the extended- 
geostrophic Hamiltonian models of A llen and Holml (Q.996 ) and A llen et al.l 
(2002}, have been derived in the Cartesian coordinates of the f3 plane 
approximation. Such approximation relies upon expansion of the equa- 
tions of motion with respect to geographic (e.g. spherical longitude and 
latitude) coordinates about some fixed point on the surface of the planet, 
in inverse powers of the (mean) radius of the planet. The expansion is 
then truncated at first order but retaining only the first order variation 
of the Coriolis parameter (the so called (3 term) and neglecting all met- 
ric terms, which are of the same order as the (3 term! Consequently, the 
(3 plane approximation is only valid locally and in geodesic coordinates 
(Phillips, 1973; Verkley, 1990). These coordinate systems are such that all 
the derivatives of the metric tensor vanishes identically at the origin and 
thus locally look like Cartesian coordinates. Geographic coordinates are 
not geodesic in general, except at the equator where coordinate curves 
are geodesic curves, e.g. great circles in spherical geometry. Consequently, 
only at the equator the (3 plane approximation is valid when written in 
geographic coordinates, but this region is forbidden for the LI model. 

Remarkable is the fact that the quasigeostrophic (QG) model — perhaps 
the most exploited model of (slow advective-timescale) intermediate dyna- 
mics — has the property of being insensitive to differences between geo- 
graphic and geodesic coordinates, namely the (3 plane approximation gives 
the right QG equations (Pedlosky, 198^1 Ripa[ I1997a| hereafter referred 
to as R97). Even though the QG system does not fit within the frame of 
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models of the LI class, i.e. it does not follow from an approximation made in 
HP for SW motion, it can be derived from HP but for stationary variations 
of a particularly chosen action ( Virasoro, 1981| Holm and Zeitlin, 1998). 



The goal of this paper is to derive an LI model using non-Cartesian 
geometry in order to make an assessment of the sensitivity of this model to 
the difference between geographic and geodesic coordinates. I am not aware 
of a similar development within the Hamiltonian framework except for the 
works of lShiittsl(|TT?89j) and |Verkley| pOOT^) . ISTTTttel (fT989|) derived a modified 



version of the Hoskins' (1975) semigeostrophic equations, which are another 
type of intermediate equations that can be derived from the LI model 



through a transformation into "geostrophic coordinates" ( Salmon, 1985 1 



Verkley (2001), in turn, presented a derivation of an isentropic Ll-type 
model for application to atmospheric flows; the model derived here is based 
on SW dynamics. Unlike both S huttst (JTp89) and Verkley s (2001J) deriva- 



tions, in this paper I use tools from non-Cartesian tensor algebra, which 
leads to an invariant formulation for the dynamical equations of the LI 
model. 

The reminder of the paper is organized as follows. In § 2 I set up a 
mathematical model for the Earth's surface that defines the space in which 
the analysis is carried out. Section 3 includes a derivation of the general 
equations for a free particle on the smooth surface of the Earth in invariant 
form. This is done from HP for a general spheroidal Earth in §3.1. The 
usual spherical approximation is then applied to the resulting motion equa- 
tions, which, in particular, are written in geographic coordinates (§3.2). 
Section 3.3 presents a discussion of the consistency of the so-called planar 
approximations, which include the classical / and (3. Section 4 is devoted 
to extending into non-Cartesian geometry Holm's (1996, hereafter referred 
to as H96) general HP for variations of Lagrangian particle labels at fixed 
Eulerian positions and time. The SW and LI model equations are derived 
in §§ 4.4 and 4.5, respectively, using the spherical Earth's model. The equa- 
tions are written in a coordinate-invariant fashion on the sphere and then 
particularized to the common geographic coordinate system. Concluding 
remarks are given in § 5. Appendix A presents various relationships involved 
in the derivation of the equations. Appendix B is reserved for the discussion 
and comparison of alternative HPs. 

2. Earth's Shape Model 

I consider here some basic geophysical facts that relate to the shape of the 
Earth and the forces acting on its equilibrium surface (e.g. lStommel and Moorel 
HlMH |Ripal 113351 |Ripa| I1997bt R97). The mathematical framework on 



which the invariant formulation of the equations derived in this paper is 
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based involve concepts from non-Cartesian tensor algebra (e.g. lAbraham et al.| 
1988: IDubrovin et al.llTM2|) that I start by reviewing first. 

2.1. NON-CARTESIAN TENSOR ALGEBRA BACKGROUND 

Let S be a two-dimensional manifold, coordinatized by x := (x 1 ,^ 2 ). Two- 
dimensional intrinsic vectors on S at any point x define the tangent space, 
TxS. The disjoint union of tangent spaces constitute the tangent bundle, 
TS. Let {ej} be a basis for T^S and {e 4 } for the dual space, (T x 5)*, namely 

e l (e 3 )=8], (2.1) 

where 5*- are the Kroenecker symbols which equal 1 if i = j and other- 
wise. Let r™(T x 5) be the space of m-contravariant and n-covariant real 
valued tensors or, simply, (m, n)-tensors. Vectors a £ Tq(T^S) = T^S are 
expressed as a = a l ei and covectors a£ T^(T^S) = (T X S)* as a = ctie 1 ; 
the quantities a 1 = a(e l ) and on = a(ei) are the components of a and a, 
respectively. (N.B. The convention of summation over repeated lower and 
upper indices is understood.) In general, a (m, n)-tensor A expresses as 

A = 4i-""i e ii ® • • • ® e im ® e jl ® • • • e jn , (2.2) 

where ...I™ = A(e 11 , • • • , e* m , ej x , • • • , e Jn ) and (g> denotes the tensor prod- 
uct. 

Assume now that S is endowed with a Riemannian metric, namely a 
symmetric, positive definite, bilinear form 

((•,•» ^roye*®^', (2.3) 

where my(x) := ((ei,ej)). The inner product of two vectors a,b £ T^S is 
computed with respect to the metric, i.e. 

((a, b)) = (rriij e 1 !g> e j ) • (a, b) = e l (a) e j (b) = rriij a % \P . (2.4) 

In particular, the square of the distance between two nearby positions on 
S, x and x + dx, is given by 

ds 2 = ((dx, dx)) = ||dx|| 2 = my dx l dx j . (2.5) 

Let b be the index lowering operator, and \ its inverse, be the index 
raising operator, which are defined by 

b : T x 5 - (T X S)*; a ^ ((a, •)) and 1 : (T X S)* - T X S, (2.6) 

respectively. The matrix of is [my], i.e. (a )j = my a J =: a«, whereas that 
of ^ is [mjj] -1 = m — 1 adj [my] =: [ra' 3 ], i.e. {a)" 1 = m^ctj =: a\ Here, 
m := det[my] and adj denotes adjoint (transpose cofactor). 
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Let, in addition, reserve the symbol d to denote the exterior derivative 
(or generalized gradient operator), whose action on a skew-symmetric (0, k)- 
tensor or fc-form a, i.e. 

a = ail ... ik e h A ■ ■ ■ A e ik , (2.7) 
where A denotes the exterior product, is defined by 

da:= d jail ... ik e j Ae h A-- - Ae ik . (2.8) 

j,h<--<i k 

Notice that, in particular, if k = then a is simply a scalar and, hence, 
da = a^e 1 =: grada. N.B. The shorthand notations di(-) and for 
partial differentiation d(-)/dx l are in use. 

Finally, let P be a linear map, with matrix elements Pjj = yjmij for i = j 
and Pjj = otherwise. Then P • a (resp., P _1 • a") denotes the physical — 
nontensorial — contravariant (resp., covariant) counterpart of vector a. For 
orthogonal coordinates, i.e. with m^- = for i / j, physical contravariant 
and covariant counterparts coincide, namely P- a = F- 1 ■ a b . 

2.2. GENERAL ASSUMPTIONS ON S 

Two main assumptions make the two-dimensional manifold S an idealized 
model of the surface of the (solid) Earth. First, S is assumed to be em- 
bedded in a three-dimensional Euclidean space which rotates steadily, with 
spinning frequency f2, with respect to a Newtonian inertial space. Second, 
S is assumed to be a geopotential surface. Namely the projections onto T X S 
of the centrifugal force (due to the spinning of the planet with respect 
to an inertial reference frame) and the gravitational attraction (due 
to the deviation of the shape of the planet from a perfect sphere and to 
inhomogeinities in the mass distribution within the planet) are assumed to 
balance one another exactly. 

As a consequence of the second assumption it follows that 

$ := V + V c = const. (2.9) 

on S, where V is the gravitational potential, Vq stands for the centrifu- 
gal potential, and their sum, $, defines the geopotential. (The constant 
in the above expression is arbitrary and can be freely set to zero.) The 
centrifugal potential (per unit mass) can be expressed in invariant form as 

^c = 4lkll 2 , (2.10) 

where a is the velocity of the x-system with respect to a suitable inertial 
frame. 
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Finally, the acceleration of gravity is denned as the minus gradient 
of <E>, thereby determining the vertical direction at each point x on S. Its 
magnitude is thus given by 



g(x) = ||grad$|| = Jm'^^j. (2.11) 



2.3. SPHERICAL MODEL 

It is convenient — and quite accurate — to consider S as a (two-dimensional) 
sphere of radius R, say, but keeping the main effect of the gravitational 
force. Namely that it can sustain a steady rotation, relative to an inertial 
frame, in any point on S. Thus let the coordinates on S be given by 

x 1 = (A - A )-Rcosi? , x 2 = (0 - i? ) R, (2-12) 

which are rescaled longitude, A, and latitude, i?, that will be referred here 
to as geographic coordinates. In this case one can introduce the usual 
notations (x, y) for (x x ,x 2 ) and (x, y) for (e 1 , e 2 ). The corresponding metric 
matrix, velocity of the x-system, and centrifugal potential, respectively, 
read: 



Here, 

/ := 2ftsintf, 7 := sec z? cos r := FT 1 tan??; (2.14) 

the first parameter is the Coriolis parameter whereas the other two are the 
geometric coefficients as defined by Ripa (2000a,b). Consistently with this 
spherical approximation, the acceleration of gravity is taken as a constant, 
namely g « 9.8 m 2 s _1 . 

More accurate models (not treated here) should account explicitly for 
the flattening of the planet at the poles. For instance, although still crude, 
next in accuracy can be mentioned one that has the form of an axisymmetric 
spheroid of revolution (Chandrasekhar, 1969; cf. also R97). 

3. Particle Dynamics 

In this section the manifold S is assumed to represent a smooth and fric- 
tionless Earth's surface on which a particle moves freely. The derivation 
of the particle's equations of motion is instructive inasmuch as it sets 
the grounds for tackling the more complicated problem of the following 
section. In particular, it shows clearly how the Coriolis force — which finds 



f - v r 



(2.13) 
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its origin in the gravitational force — arises directly from a HP with an action 
appropriate for an inertial observer, but written in coordinates fixed to the 
planet. The method is in essence the same as the one used by Pierre Simon 
de Laplace (1749-1827) to introduce this force over quarter a century before 
than Gaspard Gustave de Coriolis (1792-1843) was born (cf. R95; R96). 
The analysis of the particle's equations allows, in addition, one to simplify 
the discussion on the consistency of the so-called planar approximations 
(cf. R97). 



3.1. GENERAL EQUATIONS 



From an inertial observer viewpoint, the only force acting on the particle 
is the gravitational one. The particle's kinetic and potential energies (per 
unit mass) as measured by this observer are given by 



T(x,x) := \ ||x + <r|| 2 , y(x) = -V c = \ \W\ 



(3.1) 



respectively, where the overdot denotes time differentiation and a zero 
value of the geopotential has been assigned to the Earth's surface. The 
Lagrangian function, L : TS — ► R, is constructed in the usual way, i.e. 



L(x,x) : =T-V= i||x|| 2 + ((x,<r)). 



(3.2) 



Let 8t be a time displacement and 5x := d/de\ £=Q x(t + eSt) a variation of 
the curve x : [io,ii] — ► T X S. Let, in addition, 



S[x] :- 



dtL:F([t ,h}) 



(3.3) 



to 



be the action functional, where .FQiojii]) denotes the set of sufficiently 
smooth real valued functions on [io,ii]- Subject to fixed endpoint condi- 
tions, i.e. Sx(to) = = <$x(ti), the first variation of S, defined as 5S := 
d/de| e=0 <S[x + ei5x], is given by 



rtl / f)T 

SS= I dt ( Li5x l + — 5x' 



to 



ti 



dt L 



to 



d dL \ 



dt 



to 



\ (mjk,i ~ m ki j - m jiik ) x j x k - Xi - (a itj - a jti ) x j 5x l .(3A) 

HP (5S = 0) then yields the Newton's law for the particle in covariant form 

Dx b /dt + dcr b • x = 0. 1 (3.5) 
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In this equation, the coordinate representation of the object Da/dt, for 
any covector a, is given by (Da/di)j = dj — T^-x^ak, where r^.(x) := 

\rn kl (miij+rriji t i—mij t i) are the Christoffel symbols (of second kind), which 
establish the (Levi-Civita) connection on the Riemannian manifold S. In 
addition, 

dcr b = cTjj e l A e? = (<Jij - aj^) e l <g> e j , (3.6) 

which can be regarded as the Coriolis two-form. Notice that da ■ x = 
[(crjj - crj t i) e l <g) &>] ■ x = (dij - (Tj ti ) e'(-) e J (x) = (o it j - aj t i) x j e\ The 
operator 6 transforms (|3.5|) into its contravariant counterpart 

Dx/dt + (d<T b • x)* = 0; (3.7) 
here, Da/di, for any covector a, in components reads (Da/di) 1 = d l + 

r; 7 ,i'V\ 

Equations 1)3.5)) or 1)3.7)1 are invariant under general coordinate trans- 
formations on S, which in this case is not restricted to the spherical Earth 
model. In particular, these equations nicely show that the Coriolis term is 
responsible for the particle's trajectory to depart from a geodesic curve on 
S, i.e. a pure Galilean inertial motion. The latter is only consistent with 
motions with sufficiently large initial kinetic energy as shown by R97, who 
described all possible solutions on a sphere, namely the so-called inertial 
oscillations. 



3.2. EQUATIONS ON THE SPHERE IN GEOGRAPHIC COORDINATES 



In the geographic coordinate system ()2.12)) of the spherical Earth's model 
the only nonzero Christoffel symbols are = 7 2 r and T\ 2 = T^i = — t. 
In turn, the matrix of the Coriolis two-form takes the form 



h3 



a 



-7/ 
7/ 



(3.8) 



Thus equations 1)3.5)1 . with the spherical approximation and particularized 
to geographic coordinates, take the following component representation: 



X\ — 7 2 TX 1 X2 — (7/ + TXl)x 2 
X2 + (7/ + TXl)X l 



(3.9) 



Let u := (u, v) := P • x (= P 1 • x^ since the coordinates are orthogonal; cf. 
§2.1). Application of P _1 transforms set ()3.9)) into the more familiar form 
(e.g. R97) 

it — (/ + tu)v = 0,1 

V + (f + Tu)u = 0,J 



(3.10) 
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which can be written in vector notation as well, i.e. 

u + (/ + ru)zxu = 0, (3.11) 

where z is the vertical unit vector and x denotes the cross product of 
vectors. 

3.3. "PLANAR" APPROXIMATIONS 

In addition to the spherical approximation, other standard approximations 
introduced in the equations are the "planar" approximations. These approx- 
imations, which are meant to be valid locally at a point on the sphere in 
geographic coordinates, are obtained by expanding the equations in inverse 
powers of the radius of the sphere R. The most common approximations 
being the / and (3. The former is a consistent zeroth-order approximation. 
The latter, however, is an inconsistent first-order approximation, except at 
the equator. A consistent nth-order approximation is understood as one 
that produces 0(R^ n ^ 1 ) errors in the integrals of motion associated with 
the equations on the sphere. These integrals are the (kinetic) energy of the 
particle as measured by a terrestrial observer, 

E := ±u 2 , (3.12) 

and the absolute angular momentum (with respect to the center of the 
planet and in the direction of the axis of rotation), which, up to some 
constants, is given by 

M :=7u-OR(cost?o -7 cos??). (3.13) 

R97 showed that a consistent first-order "planar" approximation must 
have 

Ripa "plane" : 7 = 1 - r y, r = t /j, f = f + (3y/j (3.14) 

where to := -R _1 tani9o, fo ■= 20 sin -^o, and (3 := 2f2i? _1 cos #o- With 
this approximation the equations of motion conserve ^x 2 — T§yx 2 = E — 
0(R- 2 ) and (i- To y)u-f y- 1/3(1 -R 2 T$)y 2 = M-0{R- 2 ). The / plane 
approximation has 

/ plane : 7 = 1, r = 0, / = /„, (3.15) 

which consistently implies conservation of ^x 2 = E — 0(R^ 1 ) and u — foy = 
M — 0(R^ 1 ). The (3 plane approximation, in turn, has 

(3 plane : 7 = 1, r = 0, / = /„ + (3y (3.16) 
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and implies conservation of ^x 2 and u — foy — ^(3y 2 , which produce 0(ii _1 ) 
errors to E and M, respectively, everywhere except at t?q = where these 
errors are 0(R~ 2 ) because To = 0. 

It is thus clear that a consistent first-order approximation must include, 
in general, non-Cartesian terms in order to correctly reproduce the conser- 
vation laws of the system. (That is the reason for the quotation marks in 
this section.) It is worthwhile remarking that this is no longer necessary 
for motions around the equator. Geographic coordinates at the equator 
are geodesic coordinates because all the derivatives of the metric vanish 
there. For this reason locally at the equator the geometry in geographic 
coordinates looks like Cartesian and, hence, the (5 plane is a consistent 
approximation there. In general, for any point of a space with a symmet- 
ric affine connection coordinatized by x l , i = 1, 2, ••• , say, there exists 
a coordinate system x' 1 , i = 1,2, say, such that the coefficients of 
the connection vanish identically. Such a system can be defined implic- 
itly by x l = x H — ^T l j k (0) x'i x' k which can be readily seen to result in 
r^ fc (0) = 0. For geographic coordinates the transformation (x',y') ^ (x,y) 
reads (x,y) = (x',y r ) + Tox'(y', — \x'), which reduces to the identity at 
$ = i?o- Of course, the practical use of geodesic coordinates (away from the 
equator) is questionable (cf . |Phillips 119731 Verkley[ ^990). 



4. Fluid Dynamics 



In this section I derive from HP the equations of motion for (inviscid, un- 
forced) SW and LI dynamics on the spherical model for the Earth's surface. 
The derivation makes use of H96's approach but extended to non-Cartesian 
geometry. In this approach variations of Lagrangian particle labels are 
performed at fixed Eulerian positions and time. One advantage of H96's 
approach is that the equations result directly in Eulerian coordinates. 



4.1. LAGRANGIAN AND EULERIAN COORDINALES 



Identification of fluid particles in a SW motion requires two-dimensional 
labels r := (r 1 , r 2 ), say, which are defined in certain affine (metricless) space 
S, say. Let 

^xid:6xR^SxI; (j, t) i-> (x, t) = (</?(?, t),t) (4.1) 

be the map that relates the Lagrangian labels with the Eulerian two- 
dimensional positions at time t, and consider its inverse: 

ip^ 1 xidiSxl^BxR; (x, t) (r, t) = (</? -1 (x, t),t). (4.2) 
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Let now J and 3 be the Jacobians of these maps, respectively, which are 
defined by 

J := det[J?], J\ := da?/df\ (4.3a) 
3:=det[#], Z\:=df/dx\ (4.3b) 

The time derivative of a Lagrangian label, following a fluid particle, 
is zero by construction. Consequently, x = d t cp + ip^j. 1 = dfip. The latter 
defines the Lagrangian or material velocity 

t>M):=%>; (4-4) 
the Eulerian or spatial velocity, in turn, is defined by 

v(x,i):=t>(M). (4-5) 

Finally, the time derivative of any scalar function a(x, t) is a = (dt + 
v l di)a =: Da/Dt, where v 1 = v(e*). 



4.2. VOLUME CONSERVATION 



Let R(t) C S be a material spherical cap (made of the same fluid parti- 
cles) and let h(x,t) be the depth of the fluid. Let, in addition, /io(p) be 
the density of Lagrangian labels in container R(t). Since R(t) is material, 
the Lagrangian labels are defined in certain fixed region SR C 6. As a 
consequence of the metricless nature of ©, the following equality holds: 



d 2 x \/mh - 

R(t) 



dhh . 



The latter implies 



mhJ = ho, 



(4.6) 



(4.7) 



which is the Lagrangian form of the volume conservation law. In order to 
obtain the Eulerian counterpart of this law, one needs to take the time 
derivative of the l.h.s. of (|4.6j) . i.e. 



df 



d 2 x \frnh 



R(t) 



d 2 r 



d<J . — d 

Tt^ h + J <it 



mil) 



d 2 pVmJ [d t h + div(/tv)] 



(4.8) 



where the relationships of appendix have been used. The conservation 
law follows upon setting to zero the latter result: 



d t h + div(hv) = 0, 



(4.9) 
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where div(/iv) := di (yjrnhv 1 ') /y/m. Notice that in geographic coordinates 
div(hv) = 7" 1 [a r (hu) + d y (jhv)] =: V • (hu). 



4.3. GENERAL HP IN EULERIAN COORDINATES 
Following H96, I consider an action functional of the form 

£[?]:= [ 1 dtL[v,3}= ^ At [ d 2 xZ(v,5;x), (4.10) 

Jt Jt JD 

where D is a fixed region on S with solid boundary dD. Here, L is the 
Lagrangian functional and, unlike H96 who adopted Cartesian coordinates, 
l/yjm is the Lagrangian density. Variations of Lagrangian particle labels at 
fixed Eulerian positions and time result in 



5S 



D f T SL\ T SL _ , 



£5 



/ 8 < O 1 ! 



; 5L 

5v l 



(4.11) 



where /(•) := dt J £) d 2 x(-). Derivation of ()4.11j) involved the use of 
the relationships of appendix^ Fixed endpoint conditions, <5j:(x, ti) = = 
<5j(x, t2)> allows one to get rid of the last integral in (|4.11l) . Then HP implies 
the motion equation 



d t J 



r SL\ 
5v J 



+ £ v [J 



5v 



63 



and the no-flow boundary condition 



«v,n» = 



dD, 



(4.12) 



(4.13) 



where n is the external normal to the boundary. In (|4.12jh £ a a := da • a + 
d((a, a}} is the Lie derivative of covector a along vector a; in components 
(£ a a)i = a^aij + aja^. Result (|4.13j) . in turn, made use of Gauss' theorem, 
namely J D d 2 x ^Jra div a = j D d 2 x d^yjma 1 ) = §q D ds a l rii for all vector a. 
Finally, it must be mentioned that the Euler-Poincare formalism pro- 



vides an alternative way to obtaining (|4,12j) - (|4,13j) (Holm et al., 2002). 
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4.4. HP FOR SW DYNAMICS ON THE SPHERE 



Under the assumption that the layer of fluid is thin enough so that it does 
not represent a source of gravitation, an appropriate Lagrangian density 
for a HP for SW dynamics on the sphere has 



J(v,3;x) := h 3[^\W\\ 2 + ((v,a}} 



hoZ 



H 



m 



along with the definitions 



h := hoZ/ Vrn, p:= g(h- H). 



(4.14) 



(4.15) 



Here, ho and g are both constants, and p(x, t) is the hydrostatic pres- 
sure, where H(x) is the reference depth including the possibility of an 
irregular topography. The choice ho = const, is necessary in order for the 
Lagrangian density to be independent of the Lagrangian labels. The as- 
sumption g = const., in turn, is consistent with the spherical approximation 

en 



for the Earth's surface. The last term on the r.h.s. of (|4.14j) . which is not 
present in 1)3 .2|) . relates to the gravitational potential of the fluid column 
due to the departure of the free surface from the resting position. 
According to 



5L l 
7- = ho3 (v + a) , 
ov 

^ = ^(i||v|| 2 + ((v, C T))-p), 

equations (|4.12|) imply the following equivalent sets of equations: 



(4.16) 
(4.17) 



a. d t (v + a) b + £ v (v + a) b + d (p - \ ||vf - «v, a)) 



b. d t v b + d (v + a) ■ v + d [p + \ ||v| 

c. (dt + V v )v b + dcr b • v + dp = 0. 



(4.18) 



Equation (|4.18b ) involves the identity £ a a = da ■ a + d((a,a)), particu- 
larized for a = v and a = (v + a) b . Equation (14.181 ). in turn, d((a, a)) = 
VflO + V^a^- da • a — da b • oft , specialized for a = v and a = v b . Here, V a a 
denotes the covariant derivative of covector a in the direction of vector a; 
in components (V& a)i = a k ai^ — T J ik aja k . Any set selected from ()4.18jl 
together with the volume conservation equation (|4.9|1 . all subject to the 
no-flow boundary condition (|4.13|) . constitute the covariant form of the 
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SW equations on a region D denned on the sphere. These equations (or 
their contravariant counterpart via the metric) are invariant under general 
changes of coordinates on the sphere. 

The SW system conserves energy and Casimirs, namely 

S :=\J d 2 x v / ^(/i||v|| 2 + p 2 /g} , C := j d 2 x^ihC(q), (4.19) 

for arbitrary C(-) and where 



qh := -T^T-^di (j^j) = -U^ ( Vj + a,) (4.20) 



defines the potential vorticity q. The latter is conserved following fluid 
particles, i.e. T)q/Dt = 0, as readily follows upon noticing that 



*' ' dx l J-^r = (£ dx l 
dD dv l J dD 



D f T 5L\ T 5L a , 
Wt{ J 6^) +J ^ V 

dD <M 

= 0. (4.21) 

The physical counterpart of any of the equations in (|4.18j) follows from 
application of the inverse map P _1 . In geographic coordinates, the physical 
counterpart of, for instance, set (|4.18b ). reads (e.g. R97) 



d t u-(ti + f)v + 7 - 1 d x B = 0, 
dtv + (£ + f)u + d y B = 0, 



(4.22) 



where £ := , y~ l d x v — d y u — tu and B := p + \{u 2 + v 2 ) are the relative 
vorticity and Bernoulli head, respectively. System (|4.22|) can also be written 
in vector notation, i.e. 

9 iU + (( + /)zxu + VB = 0, (4.23) 

with Va := ( , y~ 1 d x a, d y a) the gradient of any scalar function a(x) in geo- 
graphic coordinates, and where £ = V • (u x z) and B = p + |u 2 . Finally, 
the integrals of motion take the form 

S = \ [ d 2 xj(hu 2 +p 2 /g) , C= [ d 2 x 7 /iC(g), (4.24) 

J D JD 

where q = (£ + f)/h satisfies 

Dq/Dt = (d t + u~f- l d x + vd y )q = (d t + u • V)q = 0. (4.25) 
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If D is a zonal channel and the topography has the same symmetry, i.e. 
d x H = 0, then zonal momentum, 

M := I d 2 x 7 /i[7u-S!i?(costfo-7costf)], (4.26) 
Jd 

is also an integral of motion. 



4.5. HP FOR LI DYNAMICS ON THE SPHERE 



The starting point of S83's method to derive approximate models by making 
approximations in the HP for SW consists is expanding the velocity field 
as 

o 7 7 + 7 <^> 

where e — > is an appropriate Rossby number. The lowest-order contribu- 
tion to the velocity is assumed to satisfy the geostrophic balance and thus 
is a function of the height (mass) field. In invariant form this reads 

v G (3) = -(d^)- 1 • dp (4.28) 

(at least there where da 9 is invertible). The Lagrangian density for LI 
dynamics on the sphere is obtained from IJ4.14J) after replacing v by (|4.27jl . 
with vq given by (|4.28j) . and by dropping the 0(e 4 )-term | ||va|| 2 in the 
first parenthesis. Thus 



Zi(v,3;x) := h 2 (((v,v G + cr)) - \ ||v G || 2 ) - \g^Jm (J^= ~ h) 



2 



(4.29) 

together with the definitions (|4.15|) . gives the LI model's Lagrangian, i.e. 
L\ := f D d 2 x l\. (A notation more consistent with my dimensional approach 
should in fact be L3 for this Lagrangian.) According to 

^ = ho3(-v G + af, (4.30) 
ov 



5L 



^- = h (((v, v G + a))-\ || v G || 2 - p AG ) , (4.31) 



LlEsfera.tex; 2/02/2008; 3:20; p. 15 



16 F. J. BERON-VERA 

where p AG '■= p — div[g , /i(d<r'') _1 • v A ]. HP implies the following equivalent 
equations: 



a. d t (v G + aY + £ v (v G + a) 9 + d (p AG - \ ||v G || 2 - ((v, v G + a)) = 0, 



2 

' (; -r lv ( ;T»] ' v -r Uj/M,; x ~ ' 

c. (dt + V v )v G + dv G • v A + dcr b • v + dp AG = 



b. 8 t v G + d (v G + af ■ v + d(p AG + i ||v G || 2 )= 0, 



(4.32) 

Because of the presence of the term ((v,v G )) in (|4.29|) . in addition to the 
no-flow boundary condition ()4.13[) . HP also implies the following condition: 



(((d^)- 1 • v A ,n)) = @ 3D. (4.33) 



Any set selected from ([4.32ft (or the corresponding contravariant counter- 
part through the metric) together with the volume conservation equation 
(|4.9j) . all subject to boundary conditions (|4,13|) and (|4.33j) . constitute the 
invariant form of the LI model on a region D on the sphere. Since v G and h 
are not independent the LI system has only one scalar prognostic equation; 
the other two scalar equations provide the constraints to determine v A . 

The LI model conserves geostrophic versions of the SW energy and 
Casimirs, namely 

£ G :=± J d 2 xv^^||v G || 2 +p 2 /5) , C G := J d 2 ^^hC(q G ) 

(4.34) 

for arbitrary C(-), where 

defines the geostrophic potential vorticity c/ G , which is materially conserved 
as is advected by the total flow (i.e. Dgc/Df = 0). 

In geographic coordinates, the physical counterpart of, for instance, set 
(|4.32b ) is given by 



dtu G - (£g + f)v + 7~ l d x B AG = 0, 
d t v G + (£g + f)u + d y B AG = 0, 



(4.36) 



where £ G := 7 - d y u G - tu g , B ag := p AG + \{v? G + v G ), and pag 

= P + 7~ 1 dx (ghv A /f) — d y (ghuA/f) — rghu A / f. In vector notation set 
(|4.36() expresses as 

d t u G + (£g + /) z x u + VB AG = 0, (4.37) 
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where ^ G = V-(u G x z) Ba g =Pag + ^u g , andp A G = p+S7-(ghu A X z/f) . 
Boundary condition (|4.33|) . in turn, takes the form 

u A • z x n = @ dD (4.38) 

(cf. |Ren and Shepherd 1 199 71 for a physical interpretation of this condition). 
The set of diagnostic equations which determines ua is given by 



A(hv A ) + B(hu A )=F 1 , 
A((g/f)hv A ) + (g/f)B(hu A ) = F 2 , 



(4.39) 



where the differential operators 

A(.) := V 2 (-) - R- 2 -r 2 - (f/g)q G , B(-) := (f / fh^d^), (4.40) 

and the functions 

Ft(h) := -d y V • (/jug) + (f/g) (hq G v G - T^xBq) , (4.41) 
F 2 {h) -^-(g/fh^dxV-ihu^ + hqGUG + dyBG; (4.42) 

here, V 2 a := V • Va = j~ 2 d xx a + d yy a — rd y a is the Laplacian of any scalar 
function a(x) in geographic coordinates*. For completeness, from (|4.39j) it 
follows 

hu A = (A(g/f)-(g/f)BA- 1 By 1 (F 2 -(g/f)BA- 1 F 1 ), (4.43a) 
hv A = (A - B(//< 7 )A- 1 ( 5 //)B)" 1 (Fx - B(f/g)A- 1 F 2 ) , (4.43b) 

which upon substitution in the volume conservation equation (|4.9[) results 
in a single evolution equation for the height field. The (Cartesian) /-plane 
version of the latter was derived by IVanneste and Bokhovel (|2UU2j) using a 
Dirac-bracket approach. Finally, the integrals of motion of the LI system 
read 

Sq := \ f d 2 x 7 {hu 2 G +p 2 /g) , C G := f d 2 x 7 hC(q G ), (4.44) 

J D JD 

where qG = (£g + f)/h; as before if D and H are zonally symmetric then 

Mg ■= / d 2 x7/i [-yu G - ^IR (cos-do ~ 7 cos??)] (4.45) 
Jd 

is also conserved. 



* Because j 2 > (excluding, of course, the poles) the elliptic problem l|4.4.'-i[l has 
a unique solution on D (bounded or periodic in one or both directions) provided that 
qg/ > ~g{R~ 2 + t 2 ) (cf. ICourant and Hilbertl IT962 ). which holds for all time because 
i/G~///»as£^0. 
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Other decompositions, appart than (|4.27jl . as well as other balance 
relationships, different than (|4.28|) . are possible (Allen and Holm, 1996; 



Alle n et al., 2002 ). This freedom is what allows for the existence of ap- 
proximate models which can be potentially more accurate than the LI 
model. 

5. Concluding Remarks 

The scaling 

{u,d t ,y/R,h-H(y),H'} = 0(e) (5.46) 
implies, at 0(e 2 ), the classical QG equation (cf. jPedlosky I1987|) 



(d t + d x ipd y - d y %l)d x ) (/qg = 0, (5.47a) 

where 

Sqg := [d xx + d yy - fi/(gH )) ^ + ((3 + fc) y. (5.47b) 

Here, H(y) = H (l- f3 T y/ f ), with H = const ., and tp := g[h-H(y)]/f is 
the geostrophic streamfunction, i.e. u = (—d y ip,d x ip) + 0(e 2 ). (More com- 
plicated topographies can of course be considered.) Notice the absence of 
geometric coefficients in (|5.47j) . Those terms, which do appear in the corre- 
sponding (diagnostic) momentum and volume conservation equations, have 
(fortuitously) cancelled out in the process of constructing the (prognostic) 
potential vorticity equation (|5.47j) [Pedlosky 1987; R97]. Consequently — 
and remarkably — QG flows develop as if the geometry were Cartesian, 
"feeling" the latitudinal variation of the Coriolis parameter as the only 
effect of the Earth's sphericity. 

The LI model shares a series of differences and similarities with the 
above QG model. Although both models are derivable from HP, the QG 
model's action is not seen to derive from approximations performed in SW's 
action. As QG motions, those governed by the LI model are not allowed at 
the equator, i.e. where / vanishes. In addition to Rossby waves, the linear 
waves of the LI model include (a form of) Kelvin waves, which are not 
supported by the QG model. Unlike QG motions, LI motions are restricted 
neither to 0{e) meridional excursions nor to 0{e) displacements of the 
free surface from the position of equilibrium at rest, nor to the presence of 
O(e) topographic variations. In a reduced-gravity setting, the equations for 
both SW and LI models have the same structure as those presented here, 
except that in that case g must be identified with the buoyancy jump at 
the interface between the active and the quiescent (infinitely deep) bottom 
layer, and H(x) must be understood as the nonuniform thickness of the 
active layer at rest, including the possibility of a nonspherical rigid surface. 
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Consequently, in contrast to the QG model, the LI model is able to describe 
the dynamics of frontal structures. 

The integrals of motion of the LI system, in geographic coordinates, 
expand in inverse powers of the radius of the spherical Earth R as 

1 _ Pv\ (d xP ? + {d y pf roy 2 _p^\ 



* = \U-% r ' h (5 - 48) 

^=((^|^-i) d -f-(^y)hy-\n^-R 2 ^)y 2 ) (s-4fl) 

C G = <(1 - roy)C(qo) + qiC(q )) (5.50) 
+ 0{R- 2 ). Here, (•) := J D d 2 x/i(-), and 

qoh:= ^±^ +f ^ (5 . 51) 
Jo 

qih := (2r ofo -P) d P -p dyyP y + r_oh_l dyp + Py . (5.52) 

Jo Jo 

Clearly, a consistent (not necessarily the optimal, though) geometric ap- 
proximation for LI dynamics, which is first-order accurate in R, is given 
by the non-Cartesian Ripa "plane" and not by the standard (3 plane (recall 
that it has to = 0). The latter is the one included in the original derivation 
of the LI model. An important contribution of the present work to the 
above list of differences and similarities between the LI and QG models 
is thus the sensitivity of the former to the differences between geographic 
and geodesic coordinates. This result confirms Ripa[s pOO Ob) in the sense 



that Earth's curvature effects increase in importance as the motions deviate 
from strictly geostrophic (divergence-free) motions. The thorough evalua- 
tion of these effects, apart from checking that the equations have the right 
conservation laws, is a subject for futher research. The latter should involve 
direct numerical simulations in which predictions of the LI model on the 
/3-plane and the sphere (or the Ripa "plane") are compared. 

Finally, |Ripa| Q1983j) showed that steady SW flows on the sphere posses 
a formal stability theorem. The latter involving an Arnold-like first theorem 
for the stability of QG flows, and a condition for the flow to be "subsonic" in 
the sense that the (geostrophic) basic flow must be everywhere slower than 
the slowest gravity- wave of the system. Ren and'Shepherd (1997) showed, 



in turn, that steady LI flows on the (3 plane posses a Ripa-like formal 
stability theorem, as well as a nonlinear (or Lyapunov) stability theorem in 
which the "subsonic" condition of Ripa's theorem is replaced by a condition 
that the flow be cyclonic along the lateral boundaries. The latter was shown 
by |Ren and Shepherdj (^997) to have an interpretation involving coastal 
Kelvin waves, which are not included in the QG model. Whether or not LI 
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flows on the sphere (or the Ripa "plane") enjoy similar stability properties 
is an issue that needs more investigation. 
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A. Useful Relations 

The following relationships can be shown to hold: 

J3 = l^ JlZ) = 5), (A.la) 

(adj J)\ = J3\^J = J*(adj J)\ <=> dJ/dJt = JZi (A. lb) 

(adja)i = ZJ\^Z = Z\{^Z)\ ^ 8Z/dZ\ = ZJ\, (A.lc) 

di (J$) = = di (24) , (A. Id) 

6 J = JdiSx\ 5Z = ZdiSf, (A.le) 

d i = j{d i ,d i = z\di. (A.if) 

In deriving IjA.lfe ) the following properties of the determinants were very 
helpful 

d(a,b) d(a,b) d(x 1 ,x 2 ) d(a,x 2 ) da 



d(t\f) d(x\x 2 ) d(t\f) ' d(x\x 2 ) 8x l 

for all scalar functions a, 6(x). 

In addition, it can be shown that: 



(A.2) 



y l = d t f + v% = => v* = -Z\d t f, (A.3a) 

dp* = Jidjj ^j{ = 4 d i v ^ ( A - 3b ) 

5v i = -J{(d t + vidjW, 5Z = ZJtdi5i\ (A.3c) 

5 J = 2d i 5x i d t Z + di (Zv 1 ) = 0. (A.3d) 
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B. Alternative HPs 

B.l. EULERIAN COORDINATES 

The standard approach for fields r(x, t) consists in considering an action 
functional of the form 

S[i]:= [ 1 dtL[ h dtt}= Pdt [ d 2 x/( f ,%,aj;x,t), (B.l) 

Jto Jto JD 

which, after invoking HP results in the familiar Euler-Lagrange equations 

(plus boundary conditions). According to di(dl/d3\) = di(3J{dl/d$) = 
3J\di(dl/dZ) and (|A.3b ). equations (|B.2|) transform into 

^( J w)- 3J ^| + ^ = °- (R3) 

The latter can be shown to be equivalent to (|4.12jl only in the particular 
case d\l = 0. 



B.2. LAGRANGIAN COORDINATES 

In the variational approach for fields x(r, t) the action is of the form 

S[x}:= [ 1 dtL[x,d t x)= [ * dt [ d 2 y /(x, d t x, J(] y, t). (B.4) 
Jt Jt is 

Upon variations of particle paths at fixed Lagrangian labels and time, the 
following Euler-Lagrange equations result from HP: 

SL 5L dl dl . . 

Ot : -jr-r = Ot—^- + O x — — — Oil = (B.5) 

5x l t ox 1 ox\ oJi 

(plus appropriate boundary conditions) . This is but the infinite-dimensional 
analogue of the particle's HP. For instance, S83's derivation of the SW and 
LI systems in the Cartesian coordinates of the /3 plane and lvan der Toornt s 
( 1997) derivation of the SW equations on the sphere are based on this HP. 
One disadvantage of this variational approach, however, is that the resulting 
equations are in Lagrangian coordinates, which requires application of the 
inverse map ip~ 1 (|4.2|) to transform back to Eulerian coordinates. 
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